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Recently the PAMELA satellite-based experiment reported an excess of galactic positrons that 
could be a signal of annihilating dark matter. The PAMELA data may admit an interpretation as 
a signal from a wino-like LSP of mass about 200 GeV, normalized to the local relic density, and 
annihilating mainly into W-bosons. This possibility requires the current conventional estimate for 
the energy loss rate of positrons be too large by roughly a factor of five. Data from anti-protons 
and gamma rays also provide tension with this interpretation, but there are significant astrophysical 
uncertainties associated with their propagation. It is not unreasonable to take this well-motivated 
candidate seriously, at present, in part because it can be tested in several ways soon. The forthcoming 
PAMELA data on higher energy positrons and the FGST (formerly GLAST) data, should provide 
important clues as to whether this scenario is correct. If correct, the wino interpretation implies a 
cosmological history in which the dark matter does not originate in thermal equilibrium. 



I. INTRODUCTION 

Recently, the PAMELA collaboration (a Pay load for 
Anti-Matter Exploration and Light-nuclei Astrophysics) 
released preliminary results [1 indicating an excess of 
cosmic ray positrons above the 10 GeV energy range. 
This confirms earlier results from HEAT [H |3] and AMS- 
01 [4j, which had already received some initial interest 
from theorists, e.g., [5l|6l[7]. 

One possible explanation for the positron excess is 
the annihilation of weakly interacting massive particles 
(WIMPS) in the galaxy. A spate of new particle physics 
models have been proposed, in part to fit the detailed fea- 
tures of the PAMELA data [SI [9]. However, it is worth 
exploring in detail whether a well-established candidate 
(such as the neutralino) could plausibly fit the data. As 
discussed recently in [T^ , for this explanation to be valid 
the neutralinos would need to have a larger cross section 
than dark matter of thermal origin. 

In light of the preliminary findings of PAMELA, we 
revisit the non-thermal neutralino models considered in 
[7J [ini [TTJ [T2] to determine whether they could account 
for the excess^. We find a wino-like neutralino with mass 
roughly 200 GeV comes close to accounting for the ex- 
cess, but only if unconventional assumptions about the 
underlying distribution of the dark matter or the prop- 
agation of its annihilation products are made. Without 
such modifications, a light supersymmetric particle ap- 
pears unable to account for the data. In part, the pur- 
pose of this paper is to point out the places where the 
sources of tension arise for this explanation, while simul- 
taneously highlighting the types of astrophysical mod- 
ifications that would need to be made to accomodate 
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the data. Such a candidate is well motivated theoreti- 
cally. For example, a wino LSP arises in theories where 
the anomaly-mediated [15 contribution to the gaugino 
masses dominates, including simple realizations of split 
supersymmetry, and the string constructions where M- 
theory is compactified on a G2 manifold |15j. 

A pure-wino 200 GeV neutralino annihilates domi- 
nantly to M^-bosons, with a cross section (av) = 2x 
10^^"* cm'^ s^^. It is remarkable that this cross section is 
approximately the correct one needed to explain the size 
of the signal in the data. Masses somewhat below 200 
GeV could conceivably explain the spectrum from the 
positron data, but such candidates come into sharper 
confiict with the existing limits from anti-protons and 
gamma rays. Even at 200 GeV, a wino has tension with 
the existing data, a fact implicit in [T6j [17]. In fact, 
taking the data at face value, such a candidate is ex- 
cluded. In the following, we will show how close the 
200 GeV neutralino comes to the current data, given the 
present understanding of the astrophysics. Given the in- 
herent astrophysical uncertainties, it is not unreasonable 
to think the 200 GeV case might ultimately be consis- 
tent with existing positron, anti-proton and 7-ray data. 
Neutralino masses much larger than this give a bad fit 
to the PAMELA results unless very large astrophysical 
boost factors are employed [IS] . 

II. POSITRON EXCESS FROM NEUTRALINO 
DARK MATTER 

Dark matter annihilations produce a differential flux 
of cosmic rays per unit time, energy and volume as 

where E is kinetic energy of the cosmic rays, {(Jxv) is the 
thermally averaged annihilation cross section and veloc- 
ity of the non-relativistic dark matter, Bi and dN^ /dE 
are the branching ratios and fragmentation functions. 
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and the sum is taken over all annihilation channels which 
contain positrons in the final state. The observed flux of 
cosmic rays is then found by propagating the source of 
Eq. (fTl), along with any astrophysical sources of cosmic 
rays (background) to the Earth. There are three sources 
of uncertainty for the prediction of any Dark Matter sig- 
nal: the dark matter distribution, the propagation of its 
annihilation products, and the role of astrophysical back- 
grounds. 

Unless otherwise stated, we restrict our attention 
to the commonly adopted Navarro-Frenk- White (NEW) 
profile |19j . This distribution is spherically symmetric 
and has the form: 

with — 20 kpc, where Tq — 8.5 kpc is the galacto- 
centric distance of the sun and pq — 0.3 GeV/cm"^ is the 
local dark matter density. Any reasonable profile will not 
effect the positrons appreciably, since they are mostly lo- 
cal. However, constraints from EGRET (and predictions 
for FGST) are directly effected by this choice, as well 
as fluxes from anti-protons. While the choice of profile 
does not appreciably effect the positron spectrum, the lo- 
cal sub-structure (dumpiness) of the dark matter, could 
have important effects, as we will soon discuss. 

For propagation of the cosmic rays resulting from 
WIMP annihilations we use DARKSUSY gD] and GAL- 
PROP [5T], to numerically solve for the fluxes. We 
vary the propagation parameters to examine how well 
a light neutralino can account for the positron excess. 
As standard values we take a diffusion coefficient oiK = 
3 X 10^'' e"'^ cm^ s~^, a half height for the conflnement re- 
gion of L = 4 kpc, and an energy loss time of r = 10^^ s. 
For the astrophysical background in positrons we adopt 
the power-law $ = (4.5e°-^) / (l GSOe^-^ -|- ISOOe^-^) 
from ,22j, where e is the energy in units of GeV. 

For comparing theoretical predictions for positrons 
with the data, it is customary to consider the positron 
fraction, which includes both contributions from dark 
matter <i>^*^ as well as the astrophysical background 

and the analogous fluxes for electrons, i.e. 

$ = sl £± (3) 

Our results for a 200 GeV particle that annihilates to 
W bosons appear in Figure [l] We have set the annihila- 
tion cross section, by assuming that the particle is a pure 
wino. The preliminary data of PAMELA are also shown 
[T], along side the anticipated signal in the presence of a 
200 GeV wino^. The bottom (solid) curve represents the 



^ Below 10 GeV, the effects of charge bias on the solar modu- 
lation are expected to be non-negligible. We have checked that 
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FIG. 1: Positron flux ratio for Wino-like Neutralino with a 
mass of 200 GeV, normalized to the local relic density. We 
set the height of the propagation region at 4 kpc and consider 
varying values for the energy loss rate (r = 1,2,5) in units 
of 10^'' s. The solid bottom line represents a conventional 
astrophysical background [22) . 

astrophysical background. The next higher curve repre- 
sents the 200 GeV wino for the NEW halo profile and de- 
fault propagation parameters discussed above. The dark 
matter signal does not provide a convincing explanation 
for the excess reported by PAMELA. 

However, we find that by varying the rate of energy 
loss of the positrons a better fit to the data is possible. 
Positrons lose energy via synchrotron radiaton and via 
inverse compoton scattering off diffuse starlight and the 
cosmic microwave background. These energy losses are 
parameterized by the energy loss time r. At the ener- 
gies of interest this is dominated by the interaction of 
the positrons with starlight. While have seen that the 
typically chosen default value of r = 10^^ s gives a poor 
fit to the data, claims in the literature |24J , indicate that 
there are theoretical uncertainties in r at the level of a 
factor of 2. So, we provide a curve for r = 2 x 10^^ s, 
still a poor fit. A t = 5 x 10^^ s gives a qualitatively 
good fit to the data. It is unclear that such a value is 
consistent with extant maps of starlight [25]. To clarify 
whether a neutralino could fit the data, it is important 
to determine this with certainty. 

The distribution of the Dark Matter could also signif- 
icantly impact the fit. Depending on the distribution of 
the dark matter (e.g., if there are significant over den- 
sities of the Dark Matter close by) astrophysical boost 
factors could also contain a dependence on the energy 



solar modulation [2JJJ can bring the PAMELA data into improved 
agreement with the theoretical estimate, but do not attempt a 
detailed accounting of the solar modulation which would require 
additional detailed data on the solar B field. As the dark mat- 
ter signal is dominant at higher energies, we focus on the region 
above 10 GeV where these effects are not important. 
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(see e.g. [25]). This would largely mimic the effects of a 
change in r, and could act to change the spectrum from 
the dark matter annihilation. An extreme example of 
this effect appeared recently in [TB] , where a local clump 
of 800 GeV wino Dark Matter was able to give the de- 
sired spectrum. We stress that the results appearing in 
Figure [T] do not include any astrophysical boost factors. 

In summary, allowing for uncertainties in the energy 
loss rate and/or allowing for a small energy dependent 
boost factor may lead to an effective value of r that could 
allow the 200 GeV candidate to account for the excess re- 
ported by PAMELA. Without invoking these uncertain- 
ties, an additional source of positrons is required. 

III. EXISTING CONSTRAINTS 

Strong bounds are set by existing data from 7-ray and 
anti- proton measurements. However both suffer from a 
number of uncertainties. 



A. Gamma Ray Constraints 

We begin with a brief review of 7-ray fluxes coming 
from dark matter annihilations, which are sensitive to 
both the halo profile and the diffuse 7-ray background. 
We then discuss existing constraints coming from the 
Energetic Gamma Ray Experiment Telescope (EGRET), 
which observed 7-rays coming from the galactic center. 

1. Overview of j-rays from dark matter Annihilation 

We are interested in the continuum energy spectrum 
of gamma-ray flux coming from neutralino annihilations. 
The differential flux is given by 

d^*^ (av) M-^dNf f o,, , , 

^ - ^E4^/i_p^(o^'W, (4) 

which is in units of photons/cm^ /s/GeV/steradian (sr). 

The first factor depends only on the particle physics. 
{av) is the thermally averaged product of the annihilation 
cross section. dNf/dE^ is the differential photon yield 
for a particular decay with branching ratio Bf, and the 
sum is taken over all relevant decays. The second piece 
contains the distribution of dark matter, where p{l) is 
the dark matter halo density profile and the integral is 
performed along the line of sight that originates from 
our location in the galaxy and continues through the full 
extent of the halo at an angle relative to the ray passing 
through the galactic-center. 

To isolate astrophysical uncertainties it is convenient 
to introduce the dimcnsionless function J{ip) 

jw = -\ I p'mm- (5) 

^QPq Jl.o.s 
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TABLE I: The averaged fine of site integral (J) to the galactic 
center for the NFW, Einasto, and Isothermal profiles with 
EGRET and FGST minimal resolution. 

Ground and satellite based detectors will observe a finite 
patch of the sky with a given angular resolution. There- 
fore, when comparing theoretical predictions with what 
may be detected, we should average J over the minimum 
angular resolution of the detector, 

where Afi is the angular resolution (in steradians). This 
value is dictated by the experiment, e.g. this corresponds 
to An = 10-3 sr for EGRET, and Aft = 10"^ sr for 
FGST. Given the minimum angular resolution, the dark 
matter profile, and the source location we can perform 
the average in (|6| using e.g. DarkSUSY. Some results 
for the line of sight integral to the galactic center appear 
in Table 1. An examination of Table 1 shows that the 
difference between a flat profile (Isothermal Cored) and 
NFW for EGRET can introduce two orders of magnitude 
difference in the signal. We also show the (J) for the 
Einasto profile, which has recently been favored by N- 
body simulations [?7] . 

While isothermal cores are now disfavored by N-body 
simulations, it is still fair to say that the current lack of 
knowledge of the halo profile induces a large error in the 
predicted fiux from the galactic center. 

Using the expression for the fiux ^ and averaging over 
the angular acceptance, the differential flux measured in 
the detector is 




which is in units of photons/cm^ /s/GeV. 

7-ray signals from dark matter annihilations must com- 
pete with the diffuse 7-ray background. These include 
inverse Compton scattering of electrons with galactic ra- 
diation and bremsstrahlung processes from accelerated 
charges |28j . Thus, uncertainties in the propagation of 
cosmic rays and in the composition of the ISM (e.g. the 
distribution and density of hydrogen) lead to uncertain- 
ties in the expected diffuse background. It is vital to 
understand the diffuse background in order to confirm 
(or deny) the existence of dark matter annihilations and 
to distinguish between different theoretical predictions. 



4 



At the present level of understanding, the differential 
flux for the diffuse 7-ray background may be fitted by a 
power-law of the form [23] 



dndK, 



(8) 



with a = -2.72 and a normalization (d2<i>^''/dridi?-y)p = 
6 X lO"'^. 



2. Constraints from EGRET 

EGRET completed nine years of observations in June 
of 2000 and was sensitive to 7-rays in the energy range 30 
MeV - 30 GeV. Using (8]) for the diffuse background near 
the galactic center and integrating over the angular res- 
olution of EGRET {An = 10"^) for the energy range of 
interest (1 GeV ^ ^ 30 GeV) we find a background 
flux of ~ 10^8 photons cm s . For dark mat- 
ter candidates that give a flux in excess of this, EGRET 
should have detected a signal. From Q, a neutralino 
annihilating to W-bosons with a mass of a couple hun- 



dred GeV and cross section (av) 



-24 



cm^ s ^ yields 



a flux comparable to the background <i> 



dm 
1 



4 X 10" 



photons cm~ s~ . This gives the first indication of the 
tension between a 200 GeV wino and the EGRET data. 
Of course, this result depends on the dark matter profile 
- assumed here to be NFW. 

Extracting robust constraints on dark matter candi- 
dates from EGRET is subtle for reasons extending be- 
yond the choice of the profile: there are uncertainties in 
both the diffuse background, as well as the EGRET data 
itself. 

EGRET has detected a possible excess above 1 GeV in 
all sky directions. Addressing the discrepancy between 
the expected diffuse background and the EGRET data 
has been considered by a number of authors [25J [301 1311 
1321 [331 [331 [3S1 [35] . These authors have argued for ex- 
planations that range from the possibility of annihilating 
dark matter^ [34] to systematic errors in the EGRET 
experiment [35]. 

A key challenge for addressing the possible excess is 
developing an accurate model of the astrophysical back- 
ground. This is particularly challenging given the in- 
ability to disentangle various components. These include 
the weak extragalactic contribution to the diffuse back- 
ground, as well as a number of possible unresolved point 
sources f3(7, '37^. Due to the uncertainties, proposed mod- 
els for the background can vary significantly. Compared 
to the background in Eqn. ([Sl, 



the 'conventional' GALPROP model [5T] assumes a 
larger contribution from inverse Compton scattering, giv- 
ing a higher contribution to the background and therefore 
to any signal that would be seen by EGRET. Yet other 
choices of background exist, including the 'Optimized' 
background [21], chosen to fit the EGRET excess with- 
out any additional Dark Matter component. At present, 
the take-home message is that there are large uncertain- 
ties in the astrophysical background. 

In addition to the uncertainties associated with the dif- 
fuse background and the halo profile, there are other rea- 
sons for concern in regards to the quality of the EGRET 
data. Indeed, EGRET was only designed to operate for 
two of its nine year mission and an aging spark chamber 
introduced time-dependent uncertainties and systematic 
errors into the high end data products [38] . In [36] is was 
found that the most likely explanation of the EGRET 
excess was an error in the estimation of the EGRET sen- 
sitivity at energies above a GeV. This was argued to be 
convincing given that the 'excess' is seen in all sky direc- 
tions, not just towards the galactic center. 

With these caveats in mind, we use the EGRET data to 
constrain the 200 GeV wino. We state the constraint as a 
bound on the (J). Assuming EGRET correctly measured 
the background above a GeV and using the data from [35] 
to determine the diffuse background, we find that a 200 
GeV wino has an annihilation cross section too large by 
a factor of three for an NFW profile - for a softer profile 
(J) ~ 380 this would not be the case. These findings 
agree with already existing bounds in the literature [301 
|37l|40]. However, we have also found that using the lower 
choice for the diffuse background in Eqn. ([S]) implies that 
the 200 GeV wino is already marginally consistent with 
the EGRET data for the NFW profile''. 

For now, it seems reasonable to consider the close prox- 
imity of the 200 GeV dark matter to the current bounds 
set by EGRET encouraging, since we will see in Section 
|IV| that the improvements of FGST should clarify the 
situation. 



B. Anti-Proton Bounds 

Measurements of cosmic ray anti-proton fluxes can also 
be used to put constraints on light neutralino candi- 
dates. In fact, the PAMELA experiment will measure 
anti-proton fluxes in the energy range 80 MeV - 190 GeV. 
It has already reported early data [42^ which seems con- 
sistent with and extends earlier results, e.g. [43l l44| |45]. 

Taken at face value, the anti-proton data would ap- 
pear to exclude a 200 GeV wino as an explanation of 
the PAMELA data, see e.g. [TBj. However, anti-proton 



^ This explanation relies on a non-standard (anisotropic) choice * See |41| for a similar approach to dealing with uncertainties in 
for the halo profile, and seems to be at odds with other sources the diffuse background and bounds on neutralinos coming from 

of indirect detection W^. EGRET and FGST. 
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FIG. 2; The anti-proton flux ratio for a 200 GeV wino- 
like neutralino as a function of kinetic energy. The lowest 
curve represents the conventional astrophysical background, 
whereas the remaining curves are the signal plus background 
for the 200 GeV candidate. These curves are the flux from 
dark matter annihilations given different choices for propaga- 
tion model - all of which have been parametrically fixed by 
matching to the well known spectrum of secondary/primary 
fluxes (e.g. B/C ratio) gg]. 



constraints suffer from theoretical uncertainties in cosmic 
ray propagation, as has been demonstrated in |46j (see 
also the discussion in [3S]). One approach to bound the 
uncertainties and set propagation parameters for anti- 
protons is to parametrically fit models to well measured 
secondary /primary fluxes such as the Boron to Carbon 
(B/C) ratio. In [iB^ it was found that this technique suf- 
fers from a number of degeneracies. These degeneracies 
arise from the fact that the effective size of the confin- 
ing region in which the cosmic rays propagate and the 
amount of energy lost to diffuse processes can be varied 
in combination, giving a good fit to the B/C ratio for a 
variety of values. 

For the 200 GeV candidate we consider here, and as- 
suming an NFW profile, the uncertainties in the propa- 
gation can lead to variations in the Dark Matter induced 
anti-proton flux by as much as an order of magnitude. 
This can be seen in Figure [2] where we present the dark 
matter signal for the benchmark propagation models ap- 
pearing in that yield the minimum and maximum 
anti-proton signal. Both models are consistent with the 
B/C ratio. The order of magnitude variation in the the- 
oretical prediction might cause the reader to be hesitant 
to conclude a 200 GeV wino is excluded from the data. 
At present, even for the minimal choice of propagation 
model, the 200 GeV candidate still gives a prediction 
that is about two times that expected from the recent 
observations of PAMELA gg. If a 200 GeV wino is to 
explain the data, there must be additional problems with 
the models used to propagate the anti-protons. 

Variation of the propagation parameters will also influ- 
ence the positron spectrum. Once the anti-proton flux is 



FIG. 3: Positron flux ratio for a wino-like neutralino with a 
mass of 200 GeV. The lowest curve represents the astrophys- 
ical background, whereas the remaining curves are the flux 
ratio for (large) energy loss rate of r = 5 x 10^® s and varying 
propagation model (as discussed in the text). 

minimized, what happens to the positrons? This effect is 
not that pronounced, primarily because the high-energy 
positrons relevant for PAMELA originate within a couple 
kpc of Earth. Propagation uncertainties are thus reduced 
relative to those for anti-protons. Fig. [3] shows the effect 
on positrons of using the "min" and "max" models used 
for Fig. [2] 

IV. FUTURE PROBES AND PREDICTIONS 
FOR FGST 

As we have mentioned, PAMELA will probe positron 
cosmic rays from 50 MeV up to an estimated 270 GeV. 
Thus, if a light wino-like neutralino is responsible for the 
positron excess, PAMELA must see a turn-over in the 
data, as can be seen from Figure [T] 

The ATIC experiment [47 has also reported an excess 
of cosmic ray electrons above roughly 300 GeV. The wino 
candidate we describe here could not be responsible for 
this excess. If the 200 GeV wino indeed accounts for the 
PAMELA excess, another explanation would be required 
for the ATIC data. 

We now consider the ability of FGST to detect the 
200 GeV wino invoked above We focus on measurements 
of the galactic center, though measurements of the halo 
could be useful if progress is made in understanding the 
backgrounds there in detail. 

FGST will offer a signiflcant improvement over 
EGRET, probing energies from 20 MeV to 300 GeV with 
an angular resolution of around 0.1 degrees (« 10^^ sr). 
The improved angular resolution will not only allow for 
separation of the point sources detected by EGRET, but 
the increased sensitivity will allow for a better opportu- 
nity to distinguish dark matter annihilation signals from 
the diffuse background. For the energy range of interest 
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FIG. 4: Tlie differential flux for the 200 GeV wino-like neu- 
tralino and an NFW profile averaged over the minimum an- 
gular resolution of FGST (i.e. AQ = 10~^ sr) and integrated 
over a 0.5° x 0.5° region around the galactic center. For the 
diffuse background we take the 'conventional' galprop model 
discussed in the text. The error bars represent the statistical 
uncertainty after one year of observations and do not account 
for systematical uncertainties. 



(around 1 GeV ^ < 300 GeV) one finds a background 
flux from ^ of around <i>Sj,^ ~ 10^^° photons cm~^ s~^ at 
a FGST angular acceptance of 10^^ sr. Compared with 
the EGRET result of ^egret ^ photons cm^^ 

s~^, this allows for an improved sensitivity by two or- 
ders of magnitude in terms of resolving signal from back- 
ground. As in the case of EGRET, the diffuse back- 
ground and halo profile are both sources of significant 
uncertainty. The better resolution and ability of FGST 
to resolve point sources should improve our knowledge of 
the diffuse background. 

For our predictions for FGST we consider an 0.5° x 0.5° 
region about the galactic center assuming an NFW profile 
and averaging with a minimum resolution set by FGST 
(i.e. 10~^ sr) . We have considered a number of choices for 
modeling the diffuse background. We find that for both 
a low choice of background given by the power-law with 
normalization in ([S]), as well as for higher backgrounds 
such as the 'conventional' and 'optimized' backgrounds 
mentioned above, that FGST will report a signal early. 
We find for the conventional background and a 200 GeV 
wino that a variation in the halo profile down to (J) ~ 70 
in the region about the galactic center can still result in a 
detectable signal for FGST at the 5cr level within the first 



year of observation. In Figure[4]we present the prediction 
for the 200 GeV wino with an NFW profile, again after 
only one year of data. The error bars reflect statistical 
uncertainties. FGST should be capable of observing the 
products of wino annihilation after the first year. 
V. CONCLUSIONS AND OUTLOOK 

The PAMELA group pLj has reported a robust excess 
in galactic positrons with energy above about 10 GeV, 
compared to those expected from known astrophysical 
sources. This is consistent with earlier reported excesses 
from [2J 131 H]. The source that produces postirons will 
likely produce antiprotons and 7-rays. No unambiguous 
excess in these channels has been observed, which pro- 
vides constraints on a dark matter interpretation. 

Only if standard assumptions about the propagation of 
cosmic rays are relaxed can the excess in positrons arise 
from galactic annihilation of winos, the superpartners of 
W bosons. Although there are apparent conflicts with 
anti-protons and gammas, given the astrophysical uncer- 
tainties it is likely premature to assume the wino inter- 
pretation is completely excluded. Further investigation 
of these uncertainties is a topic for future work. Fortu- 
nately, data from FGST and further data from PAMELA 
will help to clarify the picture soon. Later data from 
AMS-02 will also help. 

If the PAMELA excess is indeed due to a well- 
motivated wino LSP, the implications are remarkable. 
We would learn what the dark matter of our universe 
is. It would be the discovery of supersymmetry, telling 
us something about the resolution to the hierarchy prob- 
lem. It would imply a number of superpartners will likely 
be seen at LHC, confirming the result. And we would be 
learning that the universe had a non-thermal cosmologi- 
cal history that we can probe. 

Acltnowledgments 

We thank Mirko Boezio, John Beacom, Daniel Cum- 
berbatch, Jennifer Siegal-Gaskins, Igor Moskalenko, 
Brent Nelson, Malcolm Perry, Piergiorgio Picozza, Troy 
Porter, Pearl Sandick, Glenn Starkman, Greg Tarle, 
Todd Thompson, and Anna Zytkov for useful conversa- 
tions. SW would hke to thank Case Western and CCAPP 
for hospitality while this work was being completed. 

This work is supported by the Michigan Center for 
Theoretical Physics and the DOE under grant DE-FG02- 
95ER40899.The work of AP is also supported by NSF 
CAREER grant NSF-PHY-0743315. The work of SW is 
also supported by the Michigan Society of Fellows. 



[1] O. Adriani et al, arXiv:0810.4995 [astro-ph]; M. Casohno 
and P. Collaboration, arXiv:0810.4980 [astro-ph]. 

[2] S. W. Barwick et al. [HEAT Collaboration], Phys. Rev. 
Lett. 75, 390 (1995) arXiv:as tro-ph/9505141 . 

[3] S. W. Barwick et aL [HEAT Collaboration], Astrophys. 



J. 482, L191 (1997) arXiv:astro-ph/9703192'. 
[4] J. Olzem, UCLA DMS conference (2006), 

arXiv:0704.3943; A. Jacholkowska, Nucl. Phys. B 

Proc. Suppl. 165 (2007) 324. 
[5] E. A. Baltz, J. Edsjo, K. Freese and P. Gondolo, 



7 



arXiv:astro-ph/0211239 

[6] E. A. Baltz, J. Edsjo, K. Freese and P. Gondolo, Phys. 
Rev. D 65, 063511 (2002) arXiv:astro-ph/0109318 . 

[7] G. L. Kane, L. T. Wang and J. D. Wells, Phys. Rev. D 
65, 057701 (2002) arXiv:hep-ph/0108138 . 

[8] J. Hisano, M. Kawasaki, K. Kohri and K. Nakayama, 
arXiv:0812.0219 [hep-ph]; M. Ibe, H. Murayama 
and T. T. Yanagida, arXiv:0812.0072 [hep-ph]; 
M. Taoso, S. Ando, G. Bertone and S. Profumo, 
|arXiv:0811.4493| [astro-ph]; K. Ishiwata, S. Matsumoto 
and T. MoroX "arXiv:081 1.4492| [astro-ph]; K. M. Zurek, 
arXiv;0811.4429 [hep-phj; I. Cholis, G. Dobler, 

D. P. Finkbeiner, L. Goodenough and N. Weiner, 
[arXiv:08 11.3641 [astro-ph]; C. R. Chen, F. Takahashi 
and T. TT Yanagida, arXiv:0811.3357 [astro-ph]; J. Kali- 
nowski, S. F. King and J. P. Roberts, arXiv:0811.2204 
[hep-ph]; S. Baek and P. Ko, arXiv:0811.1646 [hep-ph]; 
A. Ibarra and D. Tran, arXiv:0811.1555 [hep-ph]; 
M. Pospelov, arXiv:0811.1030 [hep-ph]; E. Ponton and 
L. Randall, arXiv:0811.1029 [hep-ph]; T. Hur, H. S. Lee 
and C. Luhn, arXiv:0811.0812 [hep-ph]; K. Ham- 
aguchi, E. Nakamura, S. Shirai and T. T. Yanagida, 
|arXiv :0811.0737 [hep-ph]; C. R. Chen, F. Takahashi and 
T. T. Yanagida, arXiv:0811.0477 [hep-ph]; P. J. Fox and 

E. Poppitz, arXiv:0811.0399 [hep-ph]; Y. Bai and Z. Han, 
arXiv;0811.0387 [hep-ph]; P. f. Yin, Q. Yuan, J. Liu, 
J. Zhang, X. j. Bi and S. h. Zhu, arXiv:0811.0176 [hep- 
ph]; D. Feldman, Z. Liu and P. Nath, arXiv;0810.5762 
[hep-ph]; R. Harnik and G. D. Kribs, arXiv;0810.5557 
[hep-ph]; Y. Nomura and J. Thaler, arXiv:0810.5397 
[hep-ph]; I. Cholis, D. P. Finkbeiner, L. Goode- 
nough and N. Weiner, arXiv:0810.5344 [astro-ph]; 
T. Br ingmann, arXiv:0810.5304 [hep-ph]; P. D. Serpico, 
|arXiv :0810.4846 [hep-ph]; E. Nardi, F. Sannino and 
A. Strumia, arXiv:0811.4153 [hep-ph]. 

[9] N. Arkani-Hamed, D. P. Finkbeiner, T. Slatyer 
and N. Weiner, arXiv:0810.0713 [hep-ph]; I. Cholis, 
L. Goodenough, D. Hooper, M. Simet and N. Weiner, 
|arXiv :0809.1683 [hep-ph]; L. Bergstrom, T. Bringmann 
and J. Edsjo, arXiv:0808.3725 [astro-ph]; N. Arkani- 
Hamed and N. Weiner, arXiv:0810.0714 [hep-ph]. 
P. Grajek, G. Kane, D. J. Phalen, A. Pierce and S. Wat- 
son, arXiv:0807.1508 [hep-ph]. 

M. Nagai and K. Nakayama, Phys. Rev. D 78, 063540 
(2008) arXiv:0807.1634 [hep-ph]]'. 

T. Moroi and L. Randall, Nucl. Phys. B 570, 455 (2000) 



[20] P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schelke 
and E. A. Baltz, New Astron. Rev. 49, 149 (2005). 
P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schelke 
and E. A. Baltz, JCAP 0407, 008 (2004) |arXiv:astr(>l 



[10] 

[11] 
[12] 

[13] 
[14] 

[15] 

[16] 
[17] 
[18] 
[19] 



|arXiv :hep-ph/9906527 . 

S. 



and T. Moroi, 



K. Ishiwata, S. Matsumoto 
arXiv:0811.0250 [hep-ph]. 

L. Randall and R. Sundrum, Nucl. Phys. B 557, 
79 (1999) |arXiv:hep-th/9810155 ; G. F. Giudice, 
M. A. Luty, H. Mura yama and R. Rattazzi, JHEP 9812, 
027 (1998) |arXiv:h ep-ph/9810442|. 

B. S. Acharya, P. Kumar, K. Bobkov, G. Kane, J. Shao 
and S. Watson, JHEP 0806, 064 (2008) arXiv;0804.0863 
[hep-ph]]. 

M. CireUi, M. Kadastik, M. Raidal and A. Strumia, 
|arXiv:0809.2"409 || [hep-ph] . 

F. Donato, D. Maurin, P. Brun, T. Delahaye and 
P. Salati, arXiv:0810.5292 [astro-ph]. 
D. Hooper, A. Stebbins and K. M. Zurek, 
arXiv :0812.3202 [hep-ph]. 

J. F. Navarro, C. S. Frenk and S. D. M. White, Astro- 
phys. J. 490, 493 (1997) arXiv:astro-ph/9611107.. 



I ph/0406204 

A. W. Strong and I. V. Moskalenko, Astrophys. J. 509, 
212 (1998) |arXiv:astro-ph/9807150 . 
E. A. Baltz and J. Edsjo, Phys. Rev. D 59, 023511 (1999) 
[arXiv:astro-ph/ 9808243, . 

Clem, J. M., et. al. Astrophysical Journal v. 464, p. 507 
(1996). 

J. Olzem, H. Gast and S. Schael, Nucl. Phys. Proc. Suppl. 
173, 51 (2007). 

Private communication with Igor Moskalenko and Troy 
Porter. 

J. Lavalle, Q. Yuan, D. Maurin and X. J. Bi, 



[21 

[22 

[23; 

[24 

[25 

[26 

[27 
[28; 



[29 
[30 
[31 
[32 
[33 

[34 



[35 
[36 
[37 
[38 



arXiv:0709.3634| [astro-ph]. 

J. F. Navarro et aZ., arXiv:0810.1522 [astro-ph]. 

A. W. Strong, I. V. Moskalenko and V. S. Ptuskin, 

_Ann. Rev. Nucl. Part. Sci. 57, 285 (2007) (arXiv:astro-| 

_ph/07015l7. 

L. Bergstrom, P. Ullio and J. H. Buckley, Astropart. 
Phys. 9, 137 (1998) arXiv:astro-ph/97123l8 . 
S. Dodelson, D. Hooper and P. D. Serpico, Phys. Rev. D 
77, 063512 (2008) arXiv:071 1.4621 [astro-ph]]. 
N. Fornengo, L. Pieri and S. Scopel, Phys. Rev. D 70, 
103529 (2004) arXiv:hep-ph/0407342; . 
A. W. Strong, I. V. Moskalenko and O. Reimer, Astro- 
phys. J. 613, 962 (2004) arXiv: astro^ /0406254 . 
A. Cesarini, F. Fucito, A. Lionetto, A. Morselli and 
P. Ullio, Astropart. Phys. 21, 267 (2004) arXiv:astro^ 
ph/0305075 . ■ 
W. de Boer, C. Sander, V. Zhukov, A. V. Gladyshev and 
D. I. Kazakov, Phys. Lett. B 636, 13 (2006) |arXiv:hep-| 
_ph/0511154. 
L. Bergstrom, J. 



Edsjo, M. Gustafsson and P. Salati, 



JCAP 0605, 006 (2006) |arXiv:astro-ph/0602632 . 
F. W. Stecker, S. D. Hunter and D. A. Kniffen, Astropart. 
Phys. 29, 25 (2008) [arXlv^7d5T431 1 [astro-ph]]. 
E. A. Baltz et al., JCAP 0807, 013 (2008) 
|arXiv:0806.2911 [astro-ph]]. 

I. V. Moskalenko, S. W. Digel, T. A. Porter, O. Reimer 
and A. W. Strong, Nucl. Phys. Proc. Suppl. 173, 44 
(2007) arXiv:astro-ph/0609768 . 
[39] H. A. Mayer-Hasselwander et al, Astron. Astrophys. 

335, 161 (1998). 
[40] D. Hooper and B. L. Dingus, Phys. Rev. D 70, 113007 

(2004) arXiv;astro-ph/02106l7. 
[41] B. Altunkaynak, M. Holmes and B. D. Nelson, 

arXiv:0804.2899 [hep-ph]. 
[42] M. Boezio et al, arXiv:0810.3508 [astro-ph]; O. Adriani 

et ai, arXiv:0810.4994 [astro-ph]. 
[43] S. Orito et al. [BESS Collaboration], Phys. Rev. Lett. 84, 

1078 (2000) arXiv:astro-ph/9906426 . 
[44] T. Maeno et al. [BESS Collaboration], Astropart. Phys. 

16, 121 (2001) arXiv;astro-ph/0010381 . 
[45] K. Abe et al, arXiv:0805.1754 [astro-ph]. 
[46] F. Donato, N. Fornengo, D. Maurin and P. Salati, Phys. 

Rev. D 69, 063501 (2004) arXiv:astro-p h/0306207j . 
[47] J. Chang et al.. Nature 456, 362 (2008). 



